Observation of the antiferroquadrupolar order in DyB2C2 by resonant x-ray scattering 
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We have investigated the antiferroquadrupolar(AFQ) order in DyB2C2 by resonant x-ray scattering. 
X-rays with energies near the L3 absorption edge of Dy were employed. Superlattice peaks that 
correspond to three kinds of propagation vectors of (1 0), (1 1/2) and (0 1/2) were investigated 
in detail with polarization analyses. The experimental results are analyzed using a formalism on 
resonant x-ray scattering and a model of the AFQ order. The magnetic and quadrupolar scatterings 
are explained by the model satisfactorily. A detailed investigation on the critical behavior of the 
AFQ ordering is also reported. The critical exponent /3 is deduced to be 0.35, not far from the 
three-dimensional Heisenberg system. We have also succeeded in detecting the diffuse scattering 
above Tq. 
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I. INTRODUCTION 

There is a class of magnetic materials in which the 
orbital degree of freedom remains due to highly symmet- 
rical crystalline electric fields (CEF). Many of such ma- 
terials undergo phase transitions resulting in the lifting 
of the orbital degeneracy. In /-electron systems, when 
the CEF ground state has a non-Kramers-type degener- 
acy, there is a possibility that a periodic ordering of the 
anisotropic charge distributions of the /-electrons takes 
place, other than magnetic orderings. This asphericity 
can be represented in the lowest order by the quadrupo- 
lar moment of the localized /-electrons. The orderings 
are driven by quadrupolar pair interactions in combina- 
tion with magneto-elastic interactionsJd Even in materi- 
als in which magnetic interactions dominate, there are 
many cases where the magnetic properties are lamely 
influenced by underlying quadrupolar interactions .□ To 
study the mechanism of these phenomena the observa- 
tion of quadrupolar moments, with changing sample en- 
vironments such as temperature, pressure, and magnetic 
field is necessary. 

In order to observe the asphericity of charge distri- 
bution of an ion, diffraction method is essential, espe- 
cially for antiferroquadrupolar (AFQ) orders which nor- 
mally do not accompany measurable uniform lattice dis- 
tortions. Various methods using x-rays and neutrons 



have been studied. The most direct way to observe the 
asphericity of an electron shell is to use x-ray Thom- 
son scattering in non-resonant region. Keating's experi- 
ment on holmium succeeded in detecting the scattering 
of x-rays from the spiral arrangement of aspherical 4/- 
shells occurring at twice the magnetic propagation vec- 
tor.H The same method was formulated by Amara and 
Morin using multipole expansion, giving a description 
the scattering in terms of quadrupolar moments of ians 
The AFQ order in NdMg was observed successfully.^ We 
also refer to a x-ray powder-diffraction combined with the 
maximum entropy method which succeeded in drawing a 
charge densityf-map of a manganese oxide in an orbital 
ordered phase B 

Neutron scattering can also be used to study the 
quadrupolar orders, although neutrons do not have a 
direct coupling with electric charges. Measurement of 
magnetic form factors can give information on aspher- 
ical charge distributions through spin densities. The 
first obsesjtfation was performed by Ito and Akimitsu for 
K2CUF4.EH3 Felcher et al. measured the scattering- vector 
dependence and the temperature dependence of the as- 
pherical part of the magnetic form factor, i.e., the < j2 > 
term, of holmium by observing the third harmonic of 
the magnetic Bragg peakJj In CeBg and TmTe, the exis- 
tence of the AFQ phase with no magnetic order has been 
established through measurements of the induced anti- 
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ferromagnetic (AFM) moments in magnetic fields.ta 
Since the directions of the induced magnetic moments are 
coupled with the underlying quadrupolar moments, the 
quadrupolar order^arameter can be estimated from sym- 
metry arguments.LJ However, both the form-factor mea- 
surement and the induced-moment measurement require 
magnetic field, which changes the unperturbed state in 
zero field. 

Appropriate method for each substance and purpose 
has been employed. In holmium, since it shows ide- 
ally the spiral magnetic structure, measurement of higher 
harmonics is directly connected to the aspherical compo- 
nent. Then the detailed measurement of the temperature 
dependence of the aspherical component byjieutron scat- 
tering in zero magnetic field was possible B This is not 
the case for high symmetry systems like NdMg, CeBg 
and TmTe. We need to rely on x-ray scattering if we 
want to detect quadrupolar order in zero field, unless the 
sample exhibits periodic lattice distortion which is de- 
tectable with neutrons as in UPda.EJ A great advantage 
of the non-resonant x-ray scattering is that we can obtain 
the absolute value of the ordered quadrupolar moment. 
However, the count rate could be so weak and it might 
be difficult to collect detailed data with changing sam- 
ple environments.cHa In CeBg and TmTe, x-ray scattering 
experiments have been very difficult because of their low 
transition temperatures for x-ray experiments: Tq =3.3 
K for CeBg and 1.8 K for TmTe. Further, no indication of 
periodic lattice distortion have been observed for the two 
compounds. Then the induced AFM structures in mag- 
netic fields have been the clearest experimental evidence 
for the AFQ orders. Very recently, Nakao et al. suc- 
ceeded in observiog the AFQ order in CeBg by resonant 
x-ray scattering. Ell Finally, form-factor measurements us- 
ing x-rays and neutrons combined with the Fourier anal- 
ysis or the maximum entropy method are quite effective 
when we want to investigate actual images of spin or 
charge distributions in real space. However, these are 
not suitable for the measurement of order parameters 
with changing sample environments. 

The purpose of the present study is to observe the 
quadrupolar and magnetic order parameters of a rare- 
earth compound DyB2C2 in zero magnetic field. We have 
employed the resonant x-ray scattering for this purpose, 
which has recently been applied to the observations of 
the orbital orders in 3d transition metal oxides.Eao This 
method utilizes the characteristic that an atomic scat- 
tering factor is largely enhanced when the energy of an 
x-ray is tuned at an absorption edge of the atom. The 
first advantage of this method is the high count rate due 
to the enhancement at the absorption edge. This makes 
possible to measure temperature, azimuthal angle, polar- 
ization—and energy dependences of the peak intensities in 
detaiLo Secondly, the element selectivity guarantees the 
resonant peak arise only from the element in study with- 
out doubt. Thirdly, we can distinguish periodic arrange- 
ment of the quadrupolar moments from periodic lattice 
distortions by measuring the energy dependence of the 



superlattice peak. Finally, the experiment can be per- 
formed with a small piece of single crystal using natural 
boron, which makes neutron scattering experiment diffi- 
cult. An disadvantage is firstly that this method is not 
as direct as the non-resonant x-ray Thomson scattering; 
it is not possible to deduce the absolute value of the mo- 
ment. Secondly, the scattering mechanism has not yet 
been solidly established. Several scattering mechanisms 
in the orbital ordered state are proposed fram_different 
viewpoints, particularly in 3d compounds. Ejlj In the 
present study we simply assume that the 5d state is most 
influenced by the local 4/ electrons through the Coulomb 
and exchange interactions, which we believe is the nat- 
ural interpretation of the resonance in 4/-electron sys- 
tems. We will analyze the data quautitatively using the 
formulations developed by Blume.c3E3 

DyB2C2, with the tetragonal LaB2C2-type structure, 
is a compound that has recently been investigated in de- 
tail by Yaraauchi et al. and is considered to show an 
AFQ order .l3 Two phase transitions are clearly observed 
at 15 K and at 25 K with an entropy release of i? In 2 and 
i?ln4, respectively. The three phases are named phase 
I for T >25 K, phase II for 15 K< T <25 K, and phase 
III for T <15 K. No magnetic order appears in the phase 
II. The magnetic order in the phase HI exhibits an un- 
usual magnetic structure shown in Fig. 0. The structure 
is represented by the four propagation vectors: fci =(1 
0), k2 =(1 1/2), fes =(0 0), and =(0 1/2). The 
basic magnetic structure, where the magnetic moments 
on a c-plane lie along the [1 1 0] directions with those 
on the neighboring c-plane along the [1 1 0] directions, is 
described by fci and fc2- Canting of the moments from 
the [1 1 0] directions by an angle of 28° is described by 
fca and k^. These properties can naturally be understood 
by assuming an underlying AFQ order. The propagation 
vector of the AFQ order that is consistent with the ba- 
sic magnetic structure is expected to be (0 1/2) if we 
assume a strong spin-orbit coupling. 

The first resonant x-ray scattering |-&xperiment on 
DyB2C2 was performed by Hirota et alEB They discov- 
ered two kinds of superlattice reflections of k2 —{I 1/2) 
and fc4 =(0 1/2) that appear below Tq —25 K. From in- 
cident energy, temperature, azimuthal angle, and polar- 
ization dependences, they established that these signals 
arise from the AFQ order. In particular, the characteris- 
tic aziniuthal-angle dependences for the a— a' (oc sin^ 2ip) 
and the a — n' (cx cos^ 2ip) scatterings were considered 
to reflect the AFQ order of the 4/ electrons. Another 
reflection of fci =(1 0) was also found to appear be- 
low T/v = 15 K only for the a — n' scattering, reflecting 
the AFM order. Tanaka et al. also performed a similar 
experiment independently, though without polarizatiaa 
analysis, and obtained consistent experimental results.Eil 

This paper deals with more comprehensive data and 
analysis than the ffist reports of Refs. ^ and[3l|. The ex- 
perimental results will be connected quantitatively with 
the physical picture of the AFQ order using the basic 
formahsm described in Sec. |l[ Section HI describes the 
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experimental procedure. We will show in Sec. IV the ex- 
perimental results and the analysis in the ordered state. 
We focus especially on the explanation of the azimuthal- 
angle and the polarization dependences of the (0 5/2) 
reflection, assuming a model of the AFQ order. We also 
analyze the scatterings of magnetic origin. The critical 
phenomenon of the AFQ ordering is another topic in this 
paper, which is treated in Sec. 
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We assume here that there exists a special axis for 
the m-th atom. The axis is defined, for instance, by the 
magnetic moment, by the quadrupolar moment, or by 
the local crystalline electric field. If we take this axis as 
the a;-axis, the atomic scattering factor can be written as 



II. THEORY 

In order to analyze our experimental results we use 
the formaliS|m based on symmetry arguments developed 
by Blume.E3 We summarize the equations in this section. 
Since we deal with the scattering only in the vicinity of 
the absorpton edge, we do not consider the non-resonant 
terms. The elastic resonant scattering amplitude includ- 
ing up to electric quadrupole transition can be written 
as 
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(2) 
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are the electric dipole and quadrupole moment operators. 
fc(fc') and £(£') are the wave vector and the polarization 
vector of the incident (scattered) photon, respectively; a, 
(3, 7, and S vary over the cartesian indices x, y, and z. 
The scattering vector is written by k = k k' . The initial 
and intermediate states of the sample with energies Ea 
and Ec are represented by \a) and |c), respectively; huj is 
the energy of the photon and hcuo is equal to the energy 
difference Ec — Ea- n represents the position of the n-th 
unit cell and dm represents the m-th atom in the n-th 
unit cell with the Debye- Waller factor Wm- In Eqs. @j 
and (||) the summations are taken over all the electrons 
of the TO-th atom, pa is the thermodynamic probability 
that the sample is in the state \a). We have introduced 
r, the width of the resonance, which corresponds to the 
lifetime of the intermediate state. 

The scattering amplitude for the electric dipole (El) 
transition is written as 
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where is the atomic scattering factor tensor for the 
El transition written by 
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The parameters do, di, and c?2 are the coefficients for 
the isotropic, asymmetric, and symmetric parts of the 
tensor, respectively, including the energy dependence 
mujQ /ijj/(Tujj — huJo — ir/2). The di term arises purely 
from the magnetic moment. When the atom does not 
have a magnetic moment, the di term vanishes. The 
coefficients do and d2 contain both the magnetic and 
quadrupolar contributions. 

With respect to the scattering amplitude for the elec- 
tric quadrupole (E2) transition, we rewrite Blume's equa- 
tion in a more appealing form. We use the quadrupolar 
operators defined by 
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Then the scattering amplitude for the E2 transition is 
written as 
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where g'j^ is the atomic scattering factor tensor for the 
E2 transition written by 



{a\Qi\c){c\Q^\a) ^ 
hui — Tlloo — iT/2 



(9) 



Here the indices a and (3 vary over u, v, ^, ry, and Q. The 
coefficient K is calculated to be 
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When we assume a special axis for the m-th atom and 
take this axis as the x-axis as in the case of El transition, 
the atomic scattering factor for the E2 transition can be 
written as 



/ 9uu 


guv 















gvv 



















5?? 














V 








gvv 






























/ 








VSgvi 





\ 












gv^ 










f 






^gvi 





























9vC 






V 











-dvC 





/ 



(11) 



where 



Quu 
Quv 
Qvv 

gvv 
fee 



12&2 + 4C2 + 62 - 8/2 

-V3(4c2 + e2+4/2) 
1262 + 12c2 + 3e2 
12(62 - h) 
12(62 + C2) 
12(62 + C2) 



(12) 



and 



gv(, = -12ai 

= -12(ai + 61). 



(13) 



The coefficients 62, C2, 62, and /2 are for the symmet- 
ric part of the scattering factor and ai and 61 are for 
the asymmetric part, respectively. These are the same 
coefRcients used in Ref. 
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and include the same energy 
dependence as in the El transition. The asymmetric part 
is purely of magnetic origin and the symmetric part con- 
tains both magnetic and quadrupolar contributions. 

We notice from Eq. (|l|) that there appears the cross 
term of dipole and quadrupole transitions. We do not 
consider this contribution in this paper because the Dy 
atom is located at the center of symmetry and the cross 
term vanishes. 

The wavevectors and the polarization vectors for the 
scattering configuration illustrated in Fig. || are written 
as 



fc = A;(0, COS0, —sa\9) 
k' = fc(0, cos 6*, sin 6*) 

= (1, 0, 0) 
Ett = (0, ainO, cos 9) 

4 - (1, 0, 0) 

£^ = (0, — sin 9, cos 9). 



(14) 



When the sample is rotated around the z-axis by an 
azimuthal-angle (p, the rotation matrix 



cos (f sin ip 
U {(fi) ~ \ — sin cos If 
1 



(15) 



must be operated to all of the vectors in Eq. ( [l4| ) from 
the left. 



III. EXPERIMENTAL 

The crystal was grown by the Czochralski pulling 
method with a tetra-arc furnace. The obtained single 
crystal was checked by powder x-ray diffraction, which 
showed a single phase pattern of DyB2C2. The tem- 
perature dependence of the magnetic susceptibility also 
agreed with the data reported in Ref. 

X-ray scattering measurements were performed on 
a four circle diffractometer at BL-16A2 of the Pho- 
ton Factory in KEK. A sample with a c-plane surface 
(~2x2 mm^) was mounted in a closed cycle "^He refrig- 
erator so as to align the c-axis parallel to the (/>-axis of 
the diffractometer, i.e., the c-axis is parallel to the z-axis 
in Fig. ^. The mosaic width was 0.09° full width at half 
maximum (FWHM). The azimuthal-angle tp is defined to 
be zero when the 6-axis is parallel to the scattering plane. 
The azimuthal-angle scan can be performed by rotating 
the (/)-axis of the diffractometer. 

The incident x-ray is almost linearly polarized with its 
electric field perpendicular to the scattering plane (cr po- 
larization). The polarization of the diffracted beam, i.e., 
a' (perpendicular to the scattering plane) or tt' (parallel 
to the scattering plane), was analyzed using the PG(0 
6) reflection. The scattering angle of this reflection is 
about 91° around the L3 edge of Dy, resulting in almost 
perfect analysis. The contamination of the tt component 
in the incident x-ray was estimated to be 1.5% from the 
intensity ratio tt — tt' /a — a' of the (0 2) fundamental 
reflection. 



IV. ORDERED STATE 

A. experimental results 

The investigation of the reflection corresponding to 
fc4=(0 1/2) was performed using the (0 5/2) reflec- 
tion. The incident energy dependence of the intensity at 
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30 K, 20 K and 10 K, which correspond to the phases 
I, II, and III, respectively, is shown in Fig. |[ The mea- 
sured peak top intensities have been transformed into 
the integrated intensities by multiplying the width so as 
to be compared with other figures in the same scale. To 
compare the integrated intensities with the calculated in- 
tensities, the data have been corrected for the absorption 
and for the Lorentz factor; the absorptiosucpefhcient was 
deduced from the fluorescence spectrum£^l'L3 

Resonant peaks are clearly observed in the spectra at 
10 K and 20 K, while they are not observed at 30 K. The 
main-edge peak at 7.792 keV corresponds to the reso- 
nance due to the 2p <-> 5d electric dipole transition. The 
pre-edge peak at 7.782 keV that is well resolved in the 
a — Tr' process is probably attributed to the 2p ^ 4/ elec- 
tric quadrupole transition. These assignments are consis- 
tent with the previous experiments oa substances includ- 
ing Dy or other rare-earth elements. E3"E3 Concerning the 
a — a' process it was not possible to decide if there was 
any resonance at 7.782 keV because of the wide peak at 
the main-edge. A very small peak was observed at 10 K 
for the a — a' process at = 0°. This result is rather 
confusing since the peak is located between 7.792 keV 
and 7.782 keV. 

The temperature dependences of the integrated in- 
tensities, have been measured for each resonant peak in 
Fig. I^.Ey The resonant peaks which are observed both at 
10 K and 20 K appear below Tq = 25 K (Fig. | (b) and 
(c)), while those which are observed only at 10 K appear 
below Ttv =15 K (Fig. | (a) and (d)). Then the former 
peaks can be attributed to quadrupolar origin and the 
latter to magnetic origin. With regard to the main-edge 
and pre-edge peaks in Fig. ^ (c), the intensities at tp — 0° 
exhibit the same temperature dependence. 

Figure ^ shows the azimuthal-angle dependences of the 
integrated intensity at 7.792 keV. Peak profile was mea- 
sured by the 9 — 29 scan for each point and was fit to 
a Gaussian. The intensity of the a — a' process exhibits 
sin^ 2(/5 dependence below Tq . The intensity of the cr — tt' 
process exhibits cos^ 2lp dependence at temperatures be- 
tween Tq and T/v, and some magnetic contribution is 
added below T/v . The azimuthal-angle dependence of the 
pre-edge peak for the cr — tt' process is shown in Fig. ||. 
The intensity exhibits completely the same angle depen- 
dence with that of the main-edge, i.e., oc cos^ 2ip. With 
regard to the a — a' process, the intensity at 7.782 keV 
also exhibits the same azimuthal-angle dependence with 
that at 7.792 keV. However, it was not possible to sepa- 
rate the pre-edge component from the tail of the main- 
edge peak. 

The incident energy dependences of the integrated in- 
tensity of the (1 2) and (1 5/2) reflections, which cor- 
respond to ki and fc2, respectively, are shown in Fig. |^. 
The scale of the vertical axis is the same as those of 
Figs. ^ ^ The measurement of the temperature depen- 
dence of these resonant peaks shows that they appear be- 
low T/v, indicating magnetic origin. It should be noted 
that there is a small shoulder at the lower energy side 



of the main-edge, which is more clearly observed in the 
(1 5/2) reflection. This is considered to be of electric 
quadrupole transition. With regard to the a~a' process, 
no signal was observed at (1 2). On the other hand, the 
(1 5/2) reflection exhibits non-resmiant type energy de- 
pendence which appear below Tq.l3 This indicates that 
a periodic lattice distortion occurs simultaneously with 
the AFQ ordering. This point will be discussed in Sec. 
VI. 



B. analysis 

1. model calculation 

Let us analyze the above experimental results using 
the formalism described in Sec. ||. Figure |^ illustrates 
a model of the AFQ order of the 4/ electrons of the Dy 
ions, which is expected from the magnetic structure in 
the phase III. We introduce the canting angle a as a 
parameter. We assume each ion has its own special axis 
depending on the local quadrupolar or magnetic moment, 
around which the charge distribution is symmetric. We 
define the direction of the magnetic moment in the phase 
III as the X-axis. Since the spin-orbit coupling is strong, 
the X-axis and the principal axis of the quadrupolar mo- 
ment coincides. To calculate the scattering amplitudes 
Aei and at an azimuthal-angle ip, it is necessary to 
rotate the vectors in Eq. (^4|) so that the xyz-axes in 
Fig. 1^ coincide with those in Fig. M for each ion, i.e., the 
rotation of — 7r/4 + a + ip for Dy(l), 37r/4 — a + tp for 
Dy(2), 7r/4 + a + (p ior Dy(3), and -37r/4 - a + for 
Dy(4) is performed using Eq. (|l5|). 

The resonant scattering cross-section is equal to the 
square of the scattering amplitude A,-. The intensities of 
the (0 5/2) resonant scattering for the electric dipole 
(El) transition are calculated to be 

cx |2d2COs2asin2(/3p (16) 
l^fcr.EiP |2d2COs2acos2(^sin6' 

+2V2idi sin a sin ip cos 9\'^ , (17) 

and for the electric quadrupole (E2) transition, 

l^fe7,E2p |2(c2 + /2) cos 2a sin 2^3 sin^ 61 

+2V2{ai + bi{2 cos 2a cos 2(/3 -I- cos 2ip 
— cos2a)} cos(^sinasin20|^/16 (18) 

l^fe:,E2l' « I^{(4C2 + 62 + 4/2) + (4C2 +62+ 4/2 

-|-4c2 -I- 4/2) cos 29} cos 2a cos 2ip sin 9 
-h(terms of ai and 6i)|Vl6. (19) 

Secondly, the intensities of the (1 2) resonant scat- 
tering at = 0° for the El transition are calculated to 
be 
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(20) 
(21) 



\^k7jiif « |0.63%/2i(ii sin6'cosa|^ 
and for the E2 transition 

l^fcr.E^P |(2.7aiCos2Q! + 0.946iCos2a 

-1.26i cos 3a) sin 26*1716 (22) 
\K7^2\^ oc Kterms of ai, fei, and e2)|Vl6. (23) 

Finally, the scattering intensities of the (1 5/2) res- 
onant scattering a.t tf = 0° for the El transition are cal- 
culated to be 

\Al-'^,f^\1.9d2Sm2af (24) 
l^fcs'EiP °^ |2V2idiCosacose' + 0.5d2sin2acos6'p, (25) 
and for the E2 transition 



ia—a' |2 
^fe2,E2l 



l^fer,E2l' « Kterms of ai, &i, 62, and C2 + /2)|7l6. (27) 

The factor of |e^/mc^p has been omitted in the above 
expressions to save space. The expressions for the a — n' 
process for the E2 transition are not written explicitly 
since there appear many terms; only the parameters 
which appear are written. 



|2{0.47(c2 + /2) + 0.03162} sin 2a cos 261 
-{1.1(C2 + /2) + 0.06162} sin 2a| VI6 (26) 



2. dipole transition 

The azimuthal-angle dependences of the (0 5/2) re- 
flection at 7.792 keV in Fig. ^ are well reproduced by 
Eqs. (16) and (p7|). The sin^ 2(/5 dependence of the a — a' 
process and the cos^ 2(fi dependence of the a — tt' pro- 
cess are explained by the d2 term. Between Tjv and 
Tq, di vanishes and d2 arises only from the quadrupo- 
lar moment of the 4/-electrons. The intensity ratio 
ja-n'/jcr-a' _ g g ^^^^ explained by sin^e = 0.3155 
for the Bragg angle of the (0 5/2) reflection. The di 
term contributes below Tjy. If we assume a = 28° as de- 
termined in Ref. we obtain di = {fyz — fzy)/'2i = 0.91 
and d2 = {fxx — fyy) = 2.6 to explain quantitatively 
the integrated intensities of the (0 5/2) reflection at 
10 K by Eqs. ( p^ ) and (pT|). By comparing the inte- 
grated intensity with that of the (0 2) fundamental 
reflection (2x10'''), which is ascribed to the Thomson 
scattering from 4 x 88 electrons, we can deduce that the 
value of d2 = 2.6 corresponds to 1.4 electrons (per four 
molecules) . 

Below Tjv the ordered magnetic moment also con- 
tributes to d2 through an exchange splitting and a spin 
polarization of the 5d level.L^This should manifest below 
Tjv in the temperature dependence of the intensity. How- 
ever, the intensity for the <t — <t' process a.t ip = 45° does 



not show any clear kink in the temperature dependence 
around Tjy. This result indicates that d2 is caused mostly 
by the quadrupolar moment through Coulomb interac- 
tion between 4/ and 5d electrons. The contribution of 
the magnetic moment, which is estimated to be 7.1/is in 
Ref. seems much smaller than that of the quadrupolar 
moment. 

Using Eqs. (^) and (|25|), the above parameters at 10 
K, i.e., a = 28°, di — 0.91, and c?2 = 2.6, give the in- 
tegrated intensities of 0.11 for the (1 2) reflection and 
4.17 for the (1 5/2) reflection, respectively. The value 
of 0.11 for the (1 2) reflection does not reproduce the 
experimental result of about 1.5 in Fig. ^. This can be as- 
cribed to the experimental difficulty of rotating the crys- 
tal to different reciprocal lattice points without changing 
the effective volume which contributes to the scattering. 
However, the calculation at least explains these two mag- 
netic peaks qualitatively. 

According to Eq. (^ the intensity of the (0 5/2) 
reflection for the a ~ tt' process at ip — 45° is propor- 
tional to jdisinap. The fact that a finite intensity is 
observed at 10 K indicates that the canting angle of the 
magnetic moment is certainly not zero below Tjv. On the 
other hand, Eq. ( |25| ) shows that the intensity of the (1 
5/2) reflection for the a — tt' scattering at = 0° is 
proportional to |d2sin2ap above Tjv since the di term 
vanishes. The experimental result in Fig. |^ shows that 
the intensity completely disappears at 20 K. This indi- 
cates that the canting angle of the quadrupolar moment 
could be zero at 20 K since d2 is definitely not zero at 
this temperature. Then, it is suggested that the canting 
of the moments occurs only below T^. 

From Eq. (|2^) we see that the a — a' scattering of 
the (1 5/2) reflection can also be a measure of the 
canting angle. However, since this scattering shows a 
non-resonant type energy jiependence as a result of a 
periodic lattice distortion,E3 it was difficult to extract 
the resonant contribution at 7.792 keV with sufficient 
accuracy to examine the behavior of the canting angle in 
detail. 



3. quadrupole transition 

The calculated results are qualitatively consistent with 
the experimental results. The azimuthal-angle depen- 
dence of the (0 5/2) reflection for the a — tt' process 
at 7.782 keV shown in Fig. ^ is well reproduced by the 
quadrupolar terms of 4c2 + 62 + 4/2 = —guv/V^ and 
4(62-1-/2) = — (ffff — ff„„)/3 in Eq. (|l!^). The quadrupolar 
term of Eq. (Ilq) is also consistent with the experiment. 
The very small intensities in Fig. I (a) and (d) might be 
the magnetic signal from ai and 61 terms. 

The magnitudes of guv and {g^^ — grjri) could be esti- 
mated by comparing the intensity of the a — tt' process 
with that of the a — a' process. However, we were not 
able to estimate the pre-edge peak intensity for the a ~a' 
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process because of the difficulty in resolving the resonant 
peak into two peaks of El and E2 transitions. 

With regard to the a — it' scattering of the (1 2) and 
(1 5/2) reflections, the small shoulders around the prc- 
edge in Fig. |^ can be ascribed to the magnetic signal due 
to the oi and bi terms in Eqs. ( ^3|) and (P7|), respectively. 
The quadrupolar terms of 62 and C2 + /2 have the factor 
of sin 2a or sin4Q:. Then, the fact that the intensity 
disappears above T/v is consistent with the argument on 
the canting angle in the previous subsection. 

V. CRITICAL PHENOMENON 

Our interest here is to investigate the critical phenom- 
ena associated with the AFQ ordering by measuring the 
peak profiles precisely, especially around the transition 
temperature Tq. While almost all the orbital orderings 
in manganese oxides exhibit first-order-like phase tran- 
sitions, the AFQ ordering in DyB2C2 is considered to 
be of second order. The information on the critical phe- 
nomena may give some insight into the mechanism of the 
interaction between the quadrupole moments. 

We have utilized a Ge (111) crystal as an analyzer to 
obtain good resolution to observe diffuse scatterings. The 
longitudinal (0 I) scans for both the (0 2) fundamental 
Bragg peak and the well developed (0 5/2) superlattice 
Bragg peak at low temperature were able to be repro- 
duced by a squared Lorentzian with its FWHM 0.0011 
A~^. We have therefore assumed this as the resolution 
function for the present measurement. 

We have concentrated in this study on the (0 5/2) 
reflection a,t (p = 45°. Although the polarization analysis 
was not performed, we already know from Fig. ^that this 
reflection consists only of the a — a' scattering above T/v. 
The measurement was performed much more precisely 
than the one in Ref. |3^. Although the temperature sta- 
bility of the thermometer was kept within ±0.01 K, heat- 
ing of the sample by the beam caused a serious problem 
in the measurement of the temperature dependence of 
the weak signals at the critical region. The problem is 
that the beam intensity changes with the ring current, 
which decreases typically from 400 mA to 250 mA in a 
day. This leads to different heating powers. It reached 
about 0.25 K at most when the beam intensity was high. 
Therefore, very careful data taking and treatment of the 
data, namely shifting of the temperature, were necessary. 

The integrated intensities and the peak widths were 
obtained by fltting the profiles to a squared Lorentzian. 
Figure || shows the obtained temperature dependence of 
the integrated intensity of the (0 5/2) reflection for the 
longitudinal (0 ^) scans. Note that the scale of the verti- 
cal axis is different from the other figures. The intensities 
below Tq were fitted to a power law / cx {{Tq-T) /Tq)'^'^ 
with varying the fitting range from 15 K - 25.5 K to 23 K 
- 25.5 K. The obtained parameters are Tq = 25.52±0.009 
K and (5 — 0.35 ± 0.01, which are demonstrated by the 



solid line in the figure. This (3 is different from the pre- 
viously reported value of 0.18 in Ref. this is probably 
because the previous measurement did not have enough 
accuracy to determine the critical exponent due to small 
number of data points. 

The inset in Fig. || shows the integrated intensity 
around Tq. Although a careful data treatment was nec- 
essary to estimate the reliable temperatures as described 
above, it is certain that the intensity does not vanish 
even above Tq. Increase in the peak width was also ob- 
served around Tq . The FWHM of the original profile was 
deconvoluted with the resolution function and was con- 
verted to the inverse correlation length along the c-axis. 
The result is shown in Fig. ^. Although Kc increases with 
increasing the temperature above Tq , it was not possible 
to deduce the critical exponent from these small number 
of points. 

VI. DISCUSSIONS 
A. four propagation vectors 

We can interpret the appearance of the four propaga- 
tion vectors in the following way. At Tq = 25 K the 
AFQ order with =(0 1/2) occurs. The AFQ mo- 
ment is the principal order parameter and induces the 
periodic lattice distortion with k2 =(1 1/2) through 
some quadrupole-strain couplings. It should be noted 
that (1 1/2) and (0 1/2) are the equivalent reciprocal 
lattice points. The resonance can occur also at (1 1/2) 
simultaneously with the resonance at (0 1/2). However, 
we could not identify the resonance at (1 5/2) because 
the non-resonant scattering was dominant in the a — a' 
process or because the canting angle of the quadrupolar 
moment could be zero in the phase II. 

Below Tat = 15 K the AFM order with fci =(1 0) 
occurs. The AFM moment becomes the additional or- 
der parameter in this phase. Since (1 0) and (0 0) 
are the equivalent reciprocal lattice points, the magnetic 
Bragg peak also appears at fes =(0 0). Furthermore, 
since there is already an AFQ order which confines the 
direction of the magnetic moments, the resultant mag- 
netic structure also gives magnetic Bragg peaks at (1 
1/2) and at (0 1/2). Thus, the four propagation vectors 
are coupled with each other. 

B. periodic displacement of atoms 

The non-resonant scattering of the (1 1/2) reflection 
that appears below Tq indicates a periodic displacement 
of the atoms, which leads to the Bragg reflection due 
to the Thomson scattering. Tanaka et al. interpreted 
this reflection as a displacement of B and C atoms which, 
reduces the symmetry to the space group P42/mnm.Ell 
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Lovesey and Knight recently gave a theoretical calcula- 
tion,E3 in which they explained the (0 1/2) resonant 
peak from the same standpoint as in Ref. pl^. However, 
we consider that the displacement of B and C atoms is 
not realistic because the estimated scattering factors of 
the possible reflections from such structure are not con- 
sistent with the experimental results. If we assume the z 
parameter of the 8j site of P42 / mnm, where the B and C 
atoms are located, shifted from the original value of 0.5 to 
0.45, the squared structure factor for the (1 5/2) reflec- 
tion becomes only 0.0376 whereas those for (1 1 3/2) and 
(2 2 3/2) becomes 154.6 and 170.6, respectively. There- 
fore, much stronger reflections are expected for the two 
points. However, none of these reflections have been ob- 
served. Furthermore, neutron powder-diffraction experi- 
ments also do not show any evidence of the displacement 
of B and C atoms. CJ 

Since the atomic scattering factor of a Dy ion is much 
larger than that of B and C, it is natural to attribute 
the non-resonant (1 1/2) reflection to a periodic dis- 
placement of the Dy ions. We propose a model, where 
the Dy(l) and Dy(4) in Fig. ^ are displaced by +S along 
the c-axis and the Dy(2) and Dy(3) by —S. This gives 
the superlattice reflection at (1 1/2). The structure 
factor for the (1 5/2) reflection is calculated to be 
F = — 4j/Dy sin57r5/c, while that for the (0 2) fun- 
damental reflection is 8/b + 8/c + 4/Dy cos47r5/c. From 
the observed intensity ratio between (1 5/2) and (0 
2), a reasonable value of (5/c « 4 x 10""' is deduced. 

We consider this displacement to be related with the 
characteristic crystal structure, where the Dy ions are 
located between the hard B-C layers. This displacement 
corresponds to one of the 30 phonon modes of DyB2C2 
at k =(0 1/2). Though there are many other phonon 
modes that involve displacements of B and C, including 
the one proposed in Ref. |3^, their energy must be very 
high because they have to modify the strong covalent 
bonding among the B-C network. On the contrary, the 
energy scale of the motion of the Dy ions is expected to 
be small. Then, the position of the Dy ions is consid- 
ered to be susceptible to quadrupolar orderings through 
a coupling with the lattice. 

There is another reason that the resonant scattering in 
DyB2C2 reflects the AFQ order itself but not the sym- 
metry of the crystal. The local symmetries of Dy(l) and 
Dy(2), as numbered in Fig. ^, are different even in the 
original crystal structure above Tq. This can give the (1 
0) resonant scattering even above Tq, which is not ob- 
served experimentally. Then, even though there was any 
displacement of B and C atoms, the effect on the resonant 
scattering would be negligible. In addition, the displace- 
ment of the Dy ions as we propose is mere 4 x 10"'*, 
also negligibly small to give any effect on the resonant 
scattering. Finally, it should be noted that the theory 
in Ref. ^ can also be applied to our model of the AFQ 
order without the displacement of B and C but with the 
same space group P42 /mnm; since the theory is based 
only on the local symmetry of Dy, it is not restricted by 



the mechanism of the symmetry lowering. 

C. critical phenomenon 

It has been clarified from the present experiment that 
the order parameter vanishes continuously through Tq 
and that there is still nonvanishing intensity even above 
Tq. The temperature dependence of the order param- 
eter below Tq well follows a normal power law with a 
critical exponent (3 = 0.35 ±0.01, which is not an anoma- 
lous value when compared with the. 3D-Heisenberg model 
(0.365)Ji^ 3D-XY model (0r345),y and actual Buignetic 
systems such as EuO (0.36)^3 and MnFa (0.31).y These 
results indicate that the AFQ ordering in DyB2C2 is re- 
ally a second-order phase transition. Note that we as- 
sume in this paper that the intensity of the resonant 
scattering is proportional to the square of the order pa- 
rameter. Ishihara and Maekawa discuss this relation in 
Ref. 

Concerning the diffuse scattering, it was very diffi- 
cult to measure precisely the broadening of the width 
with increasing temperature due to the weak intensity 
in comparison with the background. However, it should 
be noted that this result does not directly mean the ab- 
sence of diffuse scattering. The present situation seems 
very similar to the case of critical pmagnetic scattering in 
MnF2 studied by x-ray scatteringJl3 It was also not pos- 
sible to measure diffuse scattering above Tn — 67.4 K 
due to weak count rates. However, the diffuse scattering 
certainly exists and were measured up to 10 K above T/v 
by neutron scattering.Cj The width of the neutron scat- 
tering profile was ten times larger than that of the x-ray 
scattering.ea 

This problem might be related with the so called two 
length scales problem. There are some systems in which 
peak profiles of diffuse scatterings above transition tem- 
peratures consist of a narrow central peak and a broad 
one. Detailed studies on this problem using both neu- 
trons and x-rays showed that x-ray signals are dominated 
by narrow components because of its small resolution vol- 
ume, ■^jl^ broad component manifests in neutron scat- 
teringoo Though the origin of the narrow component 
has not yet been established, it is interpreted as . th e near 
surface effect that is more sensitive for x-rays.E3E3 

The observed inverse correlation length Kc of DyB2C2 
above Tq is about 5 x 10~^ at a reduced tempera- 
ture (T - Tq)/Tq = 4 X 10-3, which is read from Fig. |. 
This value is as large as that of the narrow component 
of holmium measured by resonant x-ray scattering at the 
same reduced temperature (r — T/v)/lAr.c3 Furthermore, 
in holmium, broad diffuse scattering with its width ten 
times wider than the narrow component, of x-ray is cer- 
tainly observed by neutron scatteringo Therefore, we 
should not conclude that there is no broad diffuse scat- 
tering in DyB2C2 only from the present x-ray scattering 
study. 
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VII. CONCLUSIONS 

We have performed the resonant x-ray scattering on 
DyB2C2 and have investigated the signals that corre- 
spond to the AFQ and AFM orders. The experimen- 
tal results are analyzed both quantitatively and qualita- 
tively using a theory of resonant x-ray scattering and a 
model of the AFQ order. An important result is that the 
quadrupolar order parameter manifests especially in the 
a — a' scattering which little contains the magnetic con- 
tribution. We analyzed the resonant peak at (0 1/2) 
which corresponds to the AFQ order using a parameter 
d2 = {I XX - fyy) in this paper. 

The (1 1/2) reflection is also an important peak for 
more detailed study of the AFQ order in this compound. 
One reason is that the non-resonant a — a' scattering 
suggests a periodic displacement of the Dy ions below Tq, 
which is probably caused by a quadrupole-strain coupling 
that is peculiar in this compound. Another reason is 
that a detailed investigation of the resonant peak could 
reveal the behavior of the canting angle of the moments 
especially in the phase II, i.e., if it is zero or not. Both 
of those subjects require further studies. 

We have also studied the critical phenomenon of the 
AFQ ordering. The order parameter well follows a nor- 
mal power law with a reasonable critical exponent for a 
3D system. The second-order nature of the phase tran- 
sition was confirmed from the continuous decrease of the 
order parameter and the diffuse scattering above Tq. 

Finally, although the intra-atomic d — f Coulomb in- 
teraction is expected to be the most probable origin of 
the anisotropic tensor of x-ray susceptibility in the 4/- 
electron systems, we need further investigation, both the- 
oretically and experimentally, to clarify the mechanism 
of the scattering process. 
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FIG. 3. Incident energy dependences of the integrated in- 
tensity of the (0 5/2) reflections corrected for the absorption 

and the Lorentz factor: (a) o — a' scattering at <^ = 0°, (b) 
a ~ g' scattering at (p — 45° , (c) cr — tt' scattering at ip = 0° , 
and (d) a — tt' scattering aX tp = 45° . Note that the integrated 
intensity of the (0 2) fundamental peak is 2x10^. 

FIG. 4. Azimuthal-angle dependences of the integrated in- 
tensity of the (0 5/2) reflection for the a — a' and the a — ■k' 
scatterings at the main-edge. Solid lines are the flts with 
sin^ 2ip for a — a' and with cos^ 2ip for cr — tt'. 

FIG. 5. Azimuthal-angle dependence of the integrated in- 
tensity of the (0 5/2) reflection for the a — tt' scattering at 
the pre-edge. Solid line is a fit with cos^ 2(p. 

FIG. 6. Incident energy dependences of the integrated in- 
tensity for the IT — tt' scatterings a,t p = 0° corrected for the 
absorption and the Lorentz factor: (a) (10 2) reflection and 
(b) (1 5/2) reflection. 

FIG. 7. A model of the autifcrroquadrupolar order in 
DyB2C2. The shadows represent the anisotropic charge distri- 
butions. The unit cell is expressed by ax ax 2c, which contains 
four Dy ions. The canting angle a from the [1 1 0]-equivalent 
axes is treated as a parameter. The direction of the magnetic 
moment in the phase HI is taken as the x-axis. 

FIG. 8. Temperature dependence of the integrated inten- 
sity of the (0 5/2) reflection for the (0 Z)-scan at p> = 45° 
measured with a Go(lll) analyser. Solid line is a fit to a 
power law I oc {{Tq — T)/Tq)^^ . Inset shows the integrated 
intensity around the transition temperature. 

FIG. 9. Temperature dependence of the inverse correlation 
length along the c-axis obtained by deconvoluting the peak 
widths to the resolution width. Solid line is a calculated curve 
for the critical exponent f — 0.7. Inset shows the peak profile 
at T = 25.6 K. Solid line is a fit to a squared Lorentzian and 
dotted line is a squared Lorentzian with the same height and 
with the resolution width. 



FIG. 1. (left) Crystal structure of DyB2C2 {P4/mbm, 
o=5.341 A, c=3.547 A at 30 K). The magnetic structure is 
indicated by the arrows, (right) The h-l plane of the recip- 
rocal space. Black marks are the reflection points that were 
actually investigated in the present experiment. 



FIG. 2. The deflnition of the vectors associated with the 
x-rays and the axes attatched to the crystal. 
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